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In mirage-mediation models there exists a modulus field whose mass is 0(1000) TeV and its late- 
decay may significantly change the standard thermal relic scenario of the dark matter. We study 
nonthermal production of the dark matter directly from the modulus decay, and find that for some 
parameter regions non-thermally produced neutralinos can become the dark matter. 
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I. INTRODUCTION 

Supersymmetry (SUSY) [l[ is the best motivated 
physics beyond the standard model, as it provides a so- 
lution to the hierarchy problem, a natural candidate of 
the dark matter of the universe, and leads to the uni- 
fication of the gauge coupling constants. But SUSY 
must be broken spontaneously in some sector hidden 
from the observable sector, and the effect of SUSY 
breaking must be communicated to the observable sec- 
tor by some mediation mechanism. There are many 
proposed SUSY breaking and mediation mechanisms, 
such as gravity-mediation, gauge-mediation 0, H[ and 
anomaly- mediation [J]. 

On the other hand, supergravity in four space-time di- 
mensions is believed to be low-energy effective theory of 
the more fundamental higher-dimensional theory, super- 
string theory p|. One of the fundamental problems to 
obtain such a four dimensional theory from string theory 
is that there appear many moduli which are not stabi- 
lized when the extra dimensions are compactified. Also 
the light moduli cause cosmological disaster related to 
the huge entropy production from the late-time modulus 
decay pi, 0] or the runaway behavior of the dilaton || . 

Recently, Kachru-Kallosh-Lindc-Trivcdi (KKLT) Q 
provided a concrete example of the model construction 
which stabilizes all the moduli. Motivated by these re- 
cent developments, low-energy consequences of KKLT- 
typ e moduli stabilization mechanism were investigated 
[ToL HTl . [l2l [l3[ . It is found that this type of models leads 
to the mixture of modulus- and anomaly-mediation con- 
tributions to the SUSY breaking parameters, and such 
models are called mirage-mediation or mixed-modulus- 
anomaly-mediation models. Phenomenologically, there 
are some interesting features in mirage-mediation mod- 
els. One of the good points in mirage-mediation models 
is that they do not suffer from extra CP phases between 
the contributions from modulus- and anomaly-mediation. 
Another is that because of the heavy gravitino mass 
( m 3/2 ~ 0(10 — 100) TeV), the cosmological gravitino 
problem [3, [H HH may be significantly relaxed. For 
some model parameters, the little hierarchy problem can 
also be relaxed [i~7l ]. 

However, as for the cosmological issues, there remains 
other problems to be solved. In mirage-mediation mod- 



els, there exists a modulus field with mass of order 
~ (87t 2 )to3/ 2 and whose interactions with other parti- 
cles are suppressed by the Planck scale. It is likely 
that the energy stored in the coherent oscillation of the 
modulus field once dominates the universe, and reheats 
the universe with very low reheating temperature order 
~ O(100) MeV releasing huge entropy. This invalidates 
the standard thermal relic dark matter scenario, and also 
significantly dilutes the preexisting baryon asymmetry. 
Thus there is a need for investigating a consistent cosmo- 
logical scenario which can account for the present dark 
matter and baryon number of the universe in the presence 
of such a heavy modulus field. 



In order to obtain correct amount of baryon asymme- 
try, a large baryon number density should be generated 
which can even survive the dilution from the modulus de- 
cay, and perhaps the only known mechanism which can 
create such a large bar yon number is the Affleck-Dine 
mechanism [18]. In Ref. [19], it is found that Affleck-Dine 
baryogenesis mechanism can create the desired amount 
of baryon asymmetry in the presence of large dilution. 

In this paper, we study the possibility of nonthermal 
production of the dark matter directly from the modulus 
decay and find that the enhancement of the annihilation 
cross section due to the s-channel Higgs resonance makes 
it possible to account for the correct relic abundance of 
the dark matter for some parameter regions. A solu- 
tion to the moduli-induced gravitino problem [2(| is also 
discussed. A similar subject was studied in Ref. [2l[ in 
the framework of anomaly-mediated SUSY breaking sce- 
nario. Our present work can be regarded as the extension 
of their work to the mirage-mediation model, where the 
presence of moduli is necessary and inevitable. 

This paper is organized as follows. In Sec. [II] the prop- 
erties of mirage-mediation models are briefly described. 
In Sec. lIIH the cosmology with heavy moduli is discussed. 
In particular, the nonthermal dark matter from modulus 
decay and its detection possibility at on-going and future 
direct detection experiments are investigated. In Sec. IIVI 
solutions to the moduli-induced gravitino problem are 
proposed. We conclude in Sec. [Vj 
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II. A BRIEF REVIEW ON THE 
MIRAGE-MEDIATION MODEL 

In this section we briefly discuss the mirage-mediation 
models. The Kahler potential, superpotential and gauge 
kinetic fuction are assumed to be of the form (hereafter 
we set the reduced Planck scale Mp = 1) 



(H = Ai/ M , <h = mil Ml and 



K = -31n(T + T* 

W = w - Ae- aT 
fa = kT 



+ Zi(T + T*)$*$i 



(1) 



where T denotes the modulus field and denotes the 
MSSM superfields and Z ( (T + T*) = 1/(T + T*)™\ The 
subscript i distinguishes the matter species, and a and 
k are real constants. The value of rii are determined by 
the concrete model setup, but here we regard them as 
free parameters which characterize the phenomenologi- 
cal properties. This type of the Kahler potential and 
superpotential of the modulus arises from string theoret- 
ical construction such as the KKLT model. The scalar 
potential for the modulus is given by 



V s = gK(T+T* ) \ K Tf lDrWl 2 _ 3| W p 



(2) 



where D T W = W T + K T W and K TT = {K^y 1 (here 
the subscript T denotes partial derivative with it). The 
resulting minimum of this potential is a supersymmetric 
AdS vacuum. In order to obtain dS vacuum, we add the 
additional uplifting potential from the brane sequestered 
from the observable sector, 



Vuft = 



D 



(T + T*) T 



(3) 



where m is O(l) constant. Then the total scalar potential 
is given by V — V s + V\at- After fine-tuning the value 
of D, we can obtain correct dS vacuum with positive 
vacuum energy density observed now. 

This model has interesting phenomenological conse- 
quences. There are two sources of SUSY breaking in 
the present model. One comes from the uplifting brane, 
which gives the dominant contribution to the SUSY 
breaking and determines the mass of gravitino. On the 
other hand, at the dS minimum after uplifting the mod- 
ulus has non-vanishing F-term. This also gives SUSY 
breaking effect. Roughly speaking, the resulting SUSY 
breaking parameters are determined by the combination 
of anomaly-mediated SUSY breaking effect which comes 
from uplifting brane and modulus-mediated SUSY break- 
ing effect which comes from the modulus itself. 

The contribution to the gaugino masses, sfermion 
masses and A terms from the modulus-mediation at GUT 
(Grand Unified Theory) scale are denoted as M , mf and 
Ai, which are related to rij. Then we define the following 
parameters which distinguish mirage-mediation models, 



'3/2 



Mo\n(Mp/m 3/2 



(4) 



which describes the ratio of the anomaly-mediation to 
modulus-mediation contribution to the SUSY breaking 
parameters. Thus mirage-mediation models are parame- 
terized by the following parameters (T^ |. 



Mq, Cj, <Zj, tan (3, a 



(5) 



where tan/? represents the ratio of the vacuum ex- 
pectation values of the up- type and down- type Higgs, 
(H u ) I \Hd). We assume (i > 0, as positive \i is favored 
from the muon g — 2 experiment. 

When the renormalization group equations are applied 
to obtain low-energy mass spectrum, one finds that the 
gaugino masses at the energy scale /i is given by 



where 



M mW = 



g 2 a (M mir ) 



M , 



M G ut 



(M P /m 3 / 2 ) Q / 2 



(6) 



(7) 



with M GUT ~ 2 x 10 16 GeV. Thus they are unified at 
the scale fi = M m ; r , which is called the mirage-scale, 
although there are no new physics at this scale. Note that 
gauge coupling constants unification at the GUT scale 
is not modified. For a = 1, the mirage-scale actually 
becomes intermediate scale, M m ; r ~ 10 9 GeV. For a = 2, 
it becomes M m i r ~ 1 TeV. This case is called TeV scale 
mirage- mediation. 

One of the interesting properties of mirage-mediation 
models is that there exists a modulus field whose inter- 
action is Planck-suppressed. It is found that the modu- 
lus mass around its potential minimum is estimated as 
rriT ~ (87r 2 )TO 3 /2- Because the gravitino mass is one 
or two-orders of magnitude larger than the other SUSY 
particle masses, the modulus mass is expected to be 
rriT S> 0(100) TeV. The modulus with such a large mass 
decays well before Big-Bang nucleosynthesis (BBN) and 
hence it seems harmless for cosmology. But the typ- 
ical decay temperature of the modulus is much lower 
than the freeze-out temperature of the lightest supersym- 
metric particles (LSPs), which invalidates the standard 
thermal relic scenario. In the next section, we consider 
the nonthermal production of the LSPs directly from the 
modulus decay . 



1 Another problem is huge dilution of the preexisting baryon asym- 
metry due to the modulus decay. Affleck-Dine baryogenesis 
D3- E3i 23] can surely create enough baryon number which sur- 
vives the dilution Fl9l.[23. 
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III. NONTHERMAL NEUTRALINO DARK 
MATTER 

After inflation ends, the modulus is likely to dominate 
the universe before the decay because of its large energy 
density stored in the form of scalar condensates with large 
initial amplitude of order the Planck scale. Thus the final 
radiation dominated regime is realized after the modulus 
decay. The final reheating temperature is estimated as 



T T ~ 170 MeVVc ( 



rriT 



10 3 TeV 



3/2 



(8) 



where we have used the decay rate of the moduli, Ft = 
cm 3 r /47rAf|, with O(l) constant c. This satisfies the lower 
limit on the reheating temperature from BBN [25l.l26l.l27j]. 
but is well below the typical freeze-out temperature of 
the LSP, T f ~ tol/20 > O(10) GeV where m L denotes 
the mass of the LSP. Thus we can neglect the production 
process of the LSPs in the thermal bath around the epoch 
of the modulus decay. The relevant Boltzmann equations 
for the LSP, modulus and radiation are 

n L +3Hn h = -(av)nl + 2B^T T n T , (9) 
hr + 3Hut = —Tttlt, (10) 
/j rad + 4Hp lad = (m T - 2B^m h )T T n T + m^av)^ 

(11) 

where (cry) denotes the thermally averaged annihilation 
cross section of the LSP 2 , til and tit denote the number 
density of the LSP and modulus, p ra( j denotes the radia- 
tion energy density, H is the Hubble parameter, and 
denotes the branching ratio of the modulus decay into 
gauginos, which is roughly the same order as that into 
gauge bosons. For t > tt where tt is the lifetime of the 
moduli, the evolution of the LSP number density is well 
described by 



77-l + 3iJriL = — {<Jv)n L . 



(12) 



Changing the variable to Y = n-^/s where s is the en- 
tropy density, this equation can be integrated under the 
assumption that the relativistic effective degrees of free- 
dom (<7») remains constant. As a result, we obtain [30l ] 



Y{T) 



1 



Y(T T ) 



45 



(av}M P (T T -T) 



(13) 



The initial abundance of the LSP Y(Tt) is given as 

3T T 



Y(T t ) = 25* 



(14) 



2 Neutralinos are expected to soon reach kinetic equilibrium for 
Tt > O(10) MeV due to scattering with background particles in 
thermal bath [28H2I. 



Thus in the relevant parameter region, the second term in 
Eq. (|13p always dominates and the final relic abundance 
of the LSP is inversely proportional to the annihilation 
cross section similar to the usual thermal relic scenario 
[H [H, [H HE H2- The difference is that the larger 
annihilation cross section is required to obtain the desired 
amount of dark matter since the decay temperature of the 
modulus is much lower than the freeze-out temperature 
of the LSP. 

Then, whether the desired amount of dark matter is 
obtained or not crucially depends on the properties of 
the LSP. In mirage-mediation models, the lightest neu- 
tralino naturally becomes the LSP and it is the mixture 
of bino and higgsino. For small a, where the modulus- 
mediation effect dominates, the lightest neutralino is 
mainly bino-like. As the value of a is increased, the 
anomaly-mediation contribution reduces the mass of the 
gluino, which leads to lighter stop mass and smaller 
value. Thus for large a the higgsino-like LSP is also vi- 
able. 

The result is shown in Fig. [JJ for M = 350 GeV, 
cm = a M = 1, c h = o-h = and tan/3 = 10, where 
the subscript H denotes the Higgs fields and M denotes 
other matter fields. In the following analysis, we have 
used modified DarkSUSY routine [33j to calculate the an- 
nihilation cross section and spin-independent neutralino- 
nucleon scattering cross section. The purple shaded re- 
gion is excluded from the b — ► sj constraint. In the 
blue shaded region stop becomes the LSP. The gray band 
shows the favored region from the observation of WMAP 



three year results [34(. In the calculation of the relic 
abundance of the dark matter, only the S-wave contri- 
butions are included since the velocity of dark matters 
is very tiny and its effect on the annihilation cross sec- 
tion can be safely neglected. The coannihilation effects 
are also less important since the coannihilating particles 
decay very quickly compared to the rate of annihilation 
processes. Thus they are not included in our analysis. 

In general, the bino has small annihilation cross section 
and is not suitable candidate of the nonthermal dark mat- 
ter. But in mirage-mediation models there are relatively 
large parameter regions where an enhancement from the 
s-channel Higgs resonance (also called ^4-funnel) is ob- 
tained. Due to this enhancement, the bino-like neutralino 
can be the interesting candidate of the dark matter as can 
be seen from Fig. [TJ 

In Fig. [21 the spin-independent scattering cross sec- 
tion of neutralinos with nuclcons is shown for the same 
parameter as Fig. [TJ For larger a the higgsino compo- 
nent in the lightest neutralino state increases as described 
above, which leads to the enhancement on the nucleon- 
neutralino scattering cross section. The recently reported 
bound from XENON10 experiment [35[ is shown by the 
dotted line. The interesting region which accounts for the 
dark matter density of the universe satisfies the bound, 
and may be within the reach of the future direct detection 
experiments. 

More large parameter regions are investigated in Fig. [3] 
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FIG. 1: The abundance of non-thermally produced neutrali- 
nos Q. m h 2 as a function of a, for Mo = 350 GeV, cm = clm = 
1, ch = an = and tan/3 = 10. 
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FIG. 2: The spin-independent cross section of neutralinos 
with nucleons ctsi as a function of a. Parameters are the 
same as Fig. [1] 



The WMAP dark matter region is indicated by gray 
shaded region. The light gray shaded region predicts too 
low relic abundance to account for the dark matter of the 
universe, fl m h 2 < 0.1. The black upper dotted line repre- 
sents nth = 115 GeV while the lower one corresponds to 
rrih = 114 GeV, where m,h denotes the lightest Higgs bo- 



son mass. Also the neutralino-nucleon spin-independent 
scattering cross section is shown by blue dashed lines, 
csi = 10~ 45 , 10 -44 , 10~ 43 cm~ 2 from upper to lower. It 
is found that since the nonthermal scenario requires rel- 
atively light SUSY masses, the constraints from b — ► 
and Higgs mass are rather severe. In this parameter re- 
gions, the SUSY contribution for b — ► sj is dominated 
by the chargino one and it has the opposite sign to the 
standard model contribution. Since the recent theoretical 
prediction of the standard model [36[ is deviated about 
1.4 a below the experimental value 37], this chargino 
contribution severely constrains the parameters of the 
model. 
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3: The WMAP dark matter region is indicated by gray 
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0.1. The upper dotted line represents = 115 GeV while 
the lower one corresponds to = 114 GeV. Blue dashed 
lines represent asi = 10 -45 , 10 - 44 , 10 - 43 cm -2 from upper to 
lower. Other parameters are the same as Fig. [T] 



Also in Figs. [Hand [5] the results with M = 450 GeV, 
cm = clm — 1) ch = a H = and a — 1 with varying 
tan (3 are shown. Also it is found that due to the s- 
channel Higgs resonance, the resulting abundance of the 
non-thermally produced neutralinos can lie in the favored 
range. 

The contour plot is shown in Fig. [5] in the (tan j3, Mo) 
plane. In this figure a is fixed to 1. Other parame- 
ters are the same as those in Fig. |4j The WMAP dark 
matter region is indicated by gray shaded region. The 
light gray shaded region predicts too low relic abun- 
dance to account for the dark matter of the universe. 
The meanings of black dotted lines are the same as that 
of Fig. [3] Also the neutralino-nucleon spin-independent 
scattering cross section is presented by blue dashed lines, 
asi = 10~ 44 , 10~ 43 , 10~ 42 cm~ 2 from upper to lower. Al- 
though b — > 57 constraint is severe for small Mq and large 
tan /?, the s-channel resonance effect enhances the anni- 
hilation cross section which leads to smaller abundance 
of the nonthermal LSPs in rather wide parameter range. 
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FIG. 4: The abundance of non-thermally produced neu- 
tralinos Q m h 2 as a function of tan/3, for Mo = 450 GeV, 
cm = an = 1, ch = a ff = with a = 1 is fixed. 
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FIG. 6: The WMAP dark matter region is indicated by gray 
shaded region. The light gray shaded region predicts Q m h 2 < 
0.1. The upper dotted line represents nih = 115 GeV while 
the lower one corresponds to tuh = 114 GeV. Blue dashed 
lines represent asi = 10 -44 , 10 -43 , 10 -42 cm -2 from upper to 
lower. Other parameters are the same as Fig. U 
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FIG. 5: The spin-independent cross section of neutralinos 
with nucleons asi as a function of tan f3. Parameters are the 
same as Fig. U 



IV. PROBLEMS WITH NON-THERMALLY 
PRODUCED GRAVITINOS 

So far, we have ignored the possibly important pro- 
cess, the modulus decay into gravitinos. Recently it is 
pointed out that the branching ratio of the modulus de- 
cay into gravitinos -B3/2 is not suppressed and in general 



B 3/2 ~ 0(0.01 - 0.1) [Hill. The subsequent decay of 
non-thermally produced gravitinos cause another cosmo- 
logical difficulty, the destruction of light elements synthe- 
sized through BBN and overproduction of LSPs from the 
gravitino decay. Although the decay rate of moduli into 
gravitinos depends on the details of the SUSY breaking 
sector [39l . [Ifl . [4lL |4^ |. it is worth mentioning here the 
case where such a decay mode is not suppressed. 

The abundance of gravitinos produced directly from 
the modulus decay is given by 



Y, 



3/2 



2B 



3T T 



3 ^ 2 4rriT 
2.6 x 10 



9 1/2 I ^3/2 
' 0.01 



10.75 

9*(T T ) 



1/2 



V10 3 TcV/ 
(15) 



1/2 



If the gravitinos decay well before BBN, the produced 
LSPs overdose the universe. If the gravitinos decay after 
BBN starts, the decay process disturbs BBN 43]. Both 
constraints require the additional entropy production, 
which dilutes the gravitinos by an amount of A > 10 3 
after the production of gravitinos by the modulus de- 
cay. A concrete example of such an entropy produc- 
tion is the late decay of Q-balls 0, , which are non- 
topological solitons formed through Affleck-Dine mecha- 
nism [45|, Hi| [4?l [3 . This model requires no additional 
fields other than MSSM fields, and it is interesting since 
the Q-ball decay not only solves the moduli-induced grav- 
itino problem, but also can provide the correct amount 
of baryon asymmetry [l9| . 

Another possibility is thermal inflation caused by some 
additional singlet field called flaton [49|, [5(| • Thermal in- 
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flation can sufficiently dilute the modulus abundance [5l| 
in order that the subsequent modulus decay has no signif- 
icant effects on cosmology. Finally the flaton reheats the 
universe with reheating temperature of O(GeV) typically 
3 . However, it is not obvious how the sizable amount of 
baryon asymmetry is generated in the presence of ther- 
mal inflation 4 . 

After the dilution, non-thermally produced LSPs from 
the modulus decay may have neglecting abundance. 
However, in both cases, the late decay of Q-balls (or fla- 
ton) at low temperature also produces nonthermal LSPs 
[30l [551 . The resulting abundance of nonthermal LSPs is 
given by similar manner to the case of modulus decay. In 
both cases, the Q-ball (or flaton) decay into gravitinos is 
obviously suppressed since their interaction is not gravi- 
tationally suppressed, or such a decay mode can be kine- 
matically forbidden since the gravitino is heavy enough. 

In Fig. [7] the abundance of nonthermal LSPs produced 
by the Q-ball (flaton) decay for the decay temperature 
Td = 100 MeV is shown. The result is similar to that 
from the modulus decay. Note that the resulting LSP 
abundance is inversely proportional to Td, so Td can not 
be as low as 10 MeV. Contrary to the case with moduli, 
the decay temperature of the Q-ball (flaton) Td is not 
directly related to model parameters. But in order to 
dilute the gravitinos produced by the modulus decay, Td 
should not be so large. 



& m h 2 
1000 

100 
10 



T d = lOOMeV 



0.1 



0.01 



0.001 





b—>sy 


Stop 








LSP 


WMAP 

















model, motivated by the cosmological consideration that 
the modulus likely to dominate the universe before the 
decay and its decay temperature is very low. It is found 
that through the s-channel Higgs resonance, the annihila- 
tion cross section of binolike neutralino can be enhanced 
and the correct amount of dark matter can be remained. 
Although the parameter region which gives such a feature 
is constrained from the b — > sj bound and Higgs mass 
bound, the nonthermal production scenario of the dark 
matter from the modulus decay is viable for natural pa- 
rameters around a ~ 1, tan/3 ~ 10 and M Q ~ 300 — 500 
GeV. 

On the other hand, if the moduli decay into gravitinos 
with unsuppressed branching ratio, such gravitinos cause 
cosmological problems. In order to avoid the difficulty, 
the additional late-time entropy production is needed. 
Q-balls or thermal inflation flatons may be the possible 
candidates which give sufficient entropy production. We 
have also shown that although their late decay may in- 
validate the nonthermal dark matter scenario from the 
moduli, their decay itself can again produce the sizable 
amount of dark matter non-thermally. In the Q-ball sce- 
nario, not only the dark matter abundance but also the 
baryon asymmetry of the universe can be explained si- 
multaneously. 

Our final comment is that such non-thermally pro- 
duced neutralinos can not be the warm dark matter, since 
they can not freely stream without losing their energy 
due to the scattering with background par ticles. In this 
respect, the missing satellite problem [56| or the cusp 
problem [H3] should be solved without referring to the 
nature of the dark matter (see e.g., Ref. [HI). 

If any, nonthermal production scenario of the dark 
matter may have distinct implications on the collider 
and direct/indirect dark matter detection experiments 
and is interesting as a probe of the cosmological evo- 
lution scenario before BBN. Generally, in supergravity 
or superstring-inspired theory, there appear long-lived 
scalar fields whose late decay disturbs the standard ther- 
mal history of the universe. Thus we believe it is im- 
portant to seek a consistent cosmological scenario taking 
into account the dynamics of such harmful scalar fields 
in the hidden sector. 
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FIG. 7: The abundance of non-thermally produced neutrali- 



nos from the Q-ball decay Q r , 
ters are the same as Fig. [l] 



h as a function of a. Parame- 



V. CONCLUSIONS 

In this paper, we studied the possibility of nonthermal 
dark matter from the modulus decay in mirage-mediation 
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3 It is possible that the reheating temperature from the flaton de- 
cay is higher than the freeze-out temperature of the LSP. In such 
a case, the standard thermal relic scenario holds l~>2fl . 

4 See e.g., Refs. 53, 54] for related works. 
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